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Abstract 

Charge localization seems unlikely to occur in two vertically coupled symmetric quantum dots 
even if a small bias voltage breaks the exact isospin-symmetry of the system. However for a three- 
electron double quantum dot we find a strong localization of charges at certain vertically applied 
magnetic fields. The charge localization is directly connected to new ground state transitions 
between eigenstates differing only in parity. The transitions are driven by magnetic field dependent 
Coulomb correlations between the electrons and give rise to strong isospin blockade signatures in 
transport through the double dot system. 
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Quantum dot structures are excellent systems to investigate few and many particle 
physics 0, 0, [J , due to the high experimental control over the system parameters. In this 



context double quantum dots are particular 



implementation of quantum bits fqubit) j^, ls| and as a model system for molecular binding 
under controlled conditions [3, 0, Q, 0| • 

In this letter we describe a new correlation effect in a vertically coupled double quantum 
dot (DQD) in a perpendicular magnetic field, which strongly changes the molecular binding 
and at the same time, defines a two level system, that can be manipulated in a controlled 
way and could serve as a qubit. This effect is manifest in the energy spectrum and the 
transport properties of the DQD. Sweeping the magnetic field we find ground state crossings 
in a perfectly symmetric DQD containing three electrons which - in contrast to the well 
known crossings between states that differ in angular momentum and/or spin |J - occur 
between states with same spin and angular momentum. In contrast to a crossing between 
states that differ in angular momentum and/or spin, that affects the lateral motion and 
occurs in single dots already, the crossing discussed here involves a transition in the parity 
of ground state, that characterizes the vertical degree of freedom. Therefore by slightly 
breaking the symmetry between the two dots, e.g. by applying an infmitesimally small 
voltage, the crossing turns into an anticrossing, resulting in charge localization. Due to the 
charge localization, transport through the DQD is strongly suppressed at the anticrossing. 
Reducing the vertical degree of motion to an additional spin like degree of freedom, the 
isospin, the strong suppression of transport can be explained by an isospin blockade at the 
anticrossing in analogy to the well known spin blockade [lfj • 

nn 

We describe the DQD within the layer model [9|, that is applicable, if the external 
potentials separate in a strong vertical and a considerably weaker lateral component. We 
assume the in-plane confinement for the electrons to be parabolic and circular symmetric. 
Additionally a magnetic field B can be applied in the vertical direction. The in-plane motion 
of the electrons is then described in the effective mass approximation by the Fock-Darwin- 
Hamiltonian Hfd |3| and the Zeeman term Hz- 

H FD + H z ^( f+ eA)\^ + a ^ B S, (!) 

where ujq is the strength of the parabolic confinement potential, m* is the effective mass, fiB 
the Bohr Magneton and g* the effective Lande factor Q|. The eigenstates of the in-plane 



y interesting in two different manners: as an 



motion are the Fock-Darwin states \n, m) with the principal quantum number n G N and 
the angular momentum quantum number (z-component) m G Z. The Hamiltonian (JI} 
conserves the angular momentum L 2 as well as the z-component of the spin S z and the 
square of the spin S 2 , described by m, s z and s respectively. 

The vertical motion is reduced to tunneling between two 5-sheets, labeled by the quantum 
number a G {+, — }. a = ± corresponds to the upper dot (+) or lower dot (— ) respectively. 
In analogy with the real electron spin one can define a spin operator algebra, where the z- 
component of the isospin, I z is given by a [ill ]. 

The interdot tunneling Ht |±) = t |=f) which transfers electrons between the two dots 
can be expressed by isospin operators: 



H T = t + =2tl x (2) 

with the real hopping parameter i < jjj. I± are the raising and lowering operators for 
the z-component of the isospin, and I x is its x-component. The eigenstates of I x and thus 
of Ht are the symmetric and antisymmetric linear combinations of the isospin eigenstates 
|±). Due to tunneling the electrons are delocalized and the eigenstates of the Hamiltonian 
Hfd + Hp + Hz are no longer eigenstates of I z . However in case of symmetric dots the 
two layers are identical, so that the isospin-parity P is conserved. In case of more than one 
electron inside the DQD, Coulomb interaction V c between the electrons has to be included 
such that the few-electron Hamiltonian reads: 

H = H fd + Ht + Hz + V c (3) 

t=l i<j 

Since Coulomb interaction is invariant under spatial and spin rotations, total angular mo- 
mentum C z and total spin <S 2 ,S Z are still conserved and are described by the quantum 
numbers M and S, S z respectively. For a symmetric DQD also the parity of the few elec- 
tron eigenstates P = 2 Ne ■ Ix <8> • • • <8> Ix is conserved, described by the quantum number 
P G {+1,-1}. Due to Coulomb interaction the electrons are correlated. In a vertical double 
quantum dot the Coulomb interaction can be divided into two parts: Vc % ^ = V^) a + Venter- 
The intradot Coulomb interaction V^tra = ^-/ r ij describes the interaction between electrons 
localized in the same dot, whereas the interdot Coulomb interaction V^fJ r = !/(?"« + d 2 ) 1 ^ 2 
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describes the interaction between electrons localized in different dots. Here = \ri — fj\ 
is the lateral separation of two electrons % and j and d is the vertical separation between 
the dots. The Coulomb operator commutes with the z-component of the total isospin I z 
but does not commute with X x and accordingly Tij-. The commutator between Coulomb 
interaction and tunneling depends on the difference between intradot and interdot Coulomb 

n 

interaction and vanishes only in the limit d — > 0. 

Increasing the vertical magnetic field effectively leads to a stronger lateral confinement 
of the electrons and hence to an increase of the Coulomb energy. Additionally intradot 
interaction increases faster with increasing magnetic field than the interdot-interaction, that 

n. 

is limited to 1/d |9j. This different scaling causes magnetic field dependent correlations in 
the eigenstates. We show that this can lead to a ground state crossing to fixed M, S, S z for 
symmetric DQD and charge polarization in slightly asymmetric dots. 

To take correlations into account we compute the eigenstates and the corresponding 
eigenenergies by numerically diagonalizing the many-body Hamiltonian i.e. we expand 
the eigenstates in a finite basis of Slater determinants |14j . 

Calculating the magnetic-field dependence of the energy spectrum for three electrons 
inside a symmetric DQD to angular momentum M = —5 and spin S = S z = 3/2, we find 
a crossing between the two energetically lowest states as illustrated in Fig. ^ Since the 
crossing states only differ in parity, the accidental crossing converts into an anticrossing 
if the parity conservation is broken by a slight asymmetry between the dots leading to 
two strongly charge-polarized states. For specific parameters the parity crossing and hence 
the charge-polarization found for this subspace of quantum numbers becomes visible in the 
ground state (GS) as illustrated in Fig. EJ For other parameters the parity crossing will either 
arise in excited states or even dissappear completely as shown in Fig. El The asymmetry 
between the dots can be either intrinsic or caused by a small bias voltage, as it is applied 
in transport experiments Q]. We model the asymmetry between the dots by adding the 
term V z = V z ■ X z to the Hamiltonian @ where V z is the energy difference between upper 
and lower dot for a single electron. While the ground state is nearly unpolarized for general 
magnetic field strengths, Fig. El shows a strongly polarized ground state at the magnetic 
field where the anticrossing occurs. The minimal value of (2 Z ) = —0.5 corresponds to two 
electron charges in the lower dot and one in the upper. Thus we find the astonishing effect 
that electrons become localized in one of the dots by simply changing the vertical magnetic 
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FIG. 1: Energy difference of the lowest two eigenstates to 7i (M = —5 and S = S z = 3/2) as a 
function of the magnetic field B. t = —0.059 meV, Kujq = 2.96 meV and d = 19.6 nm. The crossing 
takes place at B = 7.77 T (dashed vert. line). 
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FIG. 2: Angular momentum M, total spin S and expectation value of the z-component of the 
isospin (I z ) for the three electron ground state to 1rl + V z . V z = 5.9*10 _4 X 2 meV represents a slight 
asymmetry between the dots. The peak in (I z ) illustrates the charge localization that corresponds 
to the parity crossing (see text). Other parameters are the same as in Fig. ^ 

field. It is important to note that the strength of the asymmetry (i.e. V z ) only determines 
the width of the localization dip in Fig. |2l but even for arbitrarily small asymmetries the 
ground state is strongly polarized at the anticrossing with (X z ) = —0.5. Since [C z , V z ] = 
[S 2 , V z ] = [S z , V z ] = 0, V z couples only states with same total angular momentum and total 
spin. Therefore a similar effect does not occur in the well known ground state crossing 
between states that differ in M and/or S 0. 

In the following we study the reported parity crossing in the GS of a symmetric DQD in 
more detail. Without tunneling I z is conserved and since both dots are identical the ground 
state will be two-fold degenerate with I z = ±0.5. Switching on tunneling their degeneracy 
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is lifted and the GS splits in two non-degenerate parity eigenstates \P = ±1) because of 
their different occupations of symmetric and antisymmetric orbitals. In particular Fig. ^ 
illustrates that for magnetic fields B < 7.8 T \P — —1) is favored by tunneling i.e. it has 
a higher occupation of symmetric orbitals than \P = +1). However due to magnetic-field 
dependent correlations the occupation of symmetric orbitals decreases for \P = —1) but 
increases for \P = +1), so that by increasing the magnetic field finally \P = +1) becomes 
the GS. Fig. El shows the parity as a function of tunneling and external magnetic field for 
the subspace M = —5 and spin S = S z = 3/2. The crossing exists from zero tunnel- 
ing up to t ~ 0.27 meV, which suggests to treat the tunneling t as a small perturbation. 
For small tunneling (tunneling much smaller than the energy spacing between degenerate 
ground state and first excited state at t — 0) the parity eigenstates are to first order pertur- 
bation theory given by \P = ±1) ~ (\l z = |) ± \l z = — |))/v / 2- and their energy splitting 
is 2 (l z = || 7Yt \ l z = — |) (B). As indicated this matrix element depends on the magnetic 
field due to the magnetic-field dependent correlations present in the states \l z = ±|)- To 
first order the crossing occurs at B — 7.85 T where the matrix element vanishes, and is 
independent of t in good agreement with the exact results for small tunneling (see Fig. 
For strong tunneling however higher order effects (coupling to higher states) come into play 
causing the crossing to disappear for t > 0.27 meV. Breaking the vertical symmetry of the 
DQD the two parity eigenstates are coupled and the parity-crossing converts into an anti- 
crossing, thereby lifting their accidental degeneracy by an amount V z . At the anticrossing 
the eigenstates are approximately given by I z = ±| and are thus strongly charge polar- 
ized. We want to note that the parity-crossing and the related strong charge-polarization is 
not restricted to the total angular momentum and total spin chosen here but also occurs for 
other sets of quantum numbers. Furthermore a symmetric DQD containing any odd number 
of electrons has a degenerate ground state at t = and similar parity crossings are expected 
for higher odd numbers of electrons in the DQD (the ground state of an even number of 
electrons at t = has I z = and is non-degenerate). 

The polarization of the three-electron GS can be detected in a transport experiment 
through the DQD . If a small transport voltage, Vsd, accross the DQD is ap- 

plied |l7J at constant magnetic field, the conductance G has a peak structure as a function 
of the gate voltage. The height of the conductance peaks G peak corresponding to the transi- 
tions between two and three electrons or three and four electrons inside the DQD are shown 
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FIG. 3: Dependence of the parity P for angular momentum M = —5 and S = 3/2 on magnetic 
field B and tunneling t. The solid line indicates where the crossing between the parity eigenstates 
takes place. Other parameters are the same as in Fig. ^ 

as a function of the magnetic field and for two different temperatures in Fig. 0] A com- 
parison with Fig. |21 shows that the current through the DQD is suppressed at the magnetic 
field, where the three-electron GS becomes polarized. We want to point out that since the 
asymmetry between the dots is weak only the two lowest three-electron states are polarized 
(in opposite direction) whereas the other states and in particular the two and four electron 
ground states are unpolarized (in contrast to Ref. g]). In our calculations we assume that 
transport is described by sequential tunneling processes in and out of the many-particle 
eigenstates of the isolated DQD. This is a good approximation for weak tunnel contacts 
between the reservoirs and the DQD, i.e. the tunneling strength to the external reservoirs is 
smaller than the interdot tunneling and the finite lifetime broadening of the DQD states is 
smaller than temperature . For the tunneling events a transition rate can be calculated, 
which we call T + (T~) for a transition caused by a tunneling event through the upper (lower) 
barrier. In the following we discuss the transition between two and three electrons in the 
dot, but the arguments are equally valid also for the next conductance peak. 

Assuming that an electron in the upper (lower) reservoir can only tunnel into the upper 
(lower) dot, the transition rate T + between a two particle state and a three particle state 
is proportional to the spectral weight T^ e=3 ^ Ne=2 <x J2 n ,m,a \ ( N e = 2 I d nm+a \N e = 3)| 2 
where d nm+(T denotes the annihilation operator for the orbital \nm + a) in the upper dot 
similarly =3 ^ N =2 oc J2nma K-^e = 2| d nm _ a \N e = 3)| 2 . Due to the small transport volt- 
age only the two transport channels that include the unpolarized two-electron GS and one 
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FIG. 4: Calculated height of third and fourth conductance peak (transition from N e = 2 to N e = 3 
or from N e = 3 to N e = 4)) as a function of the magnetic field for two temperatures. Parameters: 
V s d = 12^iV; T = 140 mK . Other parameters are the same as in Fig. ^ and [21 In particular the 
asymmetry between the lower and upper dot V z = 10- 2 \t\ rj 0.6/ieV. 

of the polarized three-electron states (GS and first excited state) lie within the transport 
window. Higher channels only contribute due to the finite temperature and can be further 
suppressed by lowering the temperature. For both transport channels that include polarized 
three-electron states one of the transition rates either for the tunneling in or tunneling out 
process is isospin-blocked. T~ (T + ) is suppressed, if the three electron state has two electrons 
localized in the upper (lower) dot. For a current to flow through the DQD both tunn eling 
processes are necessary, which is expressed by the effective tunneling rate T cx t-+t+ ISO!] . 
Therefore the current is strongly reduced due to an isospin blockade of both channels. Away 
from the crossing the three-electron states are no longer polarized so that the transition 
through both barriers is possible. 

We conclude: Magnetic-field dependent Coulomb correlations affect the eigenstates' tun- 
neling energies differently depending on their parity leading to additional magnetic-field 
induced level-crossings between states with different parity but same angular momentum 
and spin in a perfectly symmetric DQD. The magnetic-field dependent charge-polarization 
in symmetry-broken DQDs is due to an anticrossing of two eigenstates with different parity. 
The charge polarization takes also place in the ground state and is detectable in a transport 
experiment through the DQD as it leads to an isospin blockade at the magnetic field where 
the polarization occurs. The resulting polarized eigenstates ||) and \ — \) can be seen as a 
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qubit which can be switched by the applied bias voltage. A controlled superposition of the 
two states can then be achieved by adjusting the magnetic field. The parity crossing and 
the related magnetic-field dependent charge polarization appear also for an odd number of 
electrons greater than three and in different subsets of quantum numbers. 

We thank Michael Tews for his intensive help with the transport calculations and ac- 
knowledge financial support by the Deutsche Forschungsgemeinschaft via SFB 508. 
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